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Abstract 



Searches for resonant v production in e + e~ collisions under the assumption that 
R-parity is not conserved and that the dominant R-parity violating coupling is 
A121 or A131 used data recorded by DELPHI in 1997 to 2000 at centre-of-mass 
energies of 183 to 208 GeV. No deviation from the Standard Model was ob- 
served. Upper limits are given for the Ai 2 i and Ai 3 i couplings as a function of 
the sneutrino mass and total width. The limits are especially stringent for sneu- 
trino masses equal to the centre-of-mass energies with the highest integrated 
luminosities recorded. 
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1 Introduction 

In the Minimal Supersymmetric extension of the Standard Model (MSSM) [1], a dis- 
crete symmetry called R-parity [2] predicts baryon number (B) and lepton number (L) 
conservation which is an accidental feature of the SU(3)xSU(2)xU(l) Standard Model 
(SM). The related quantum number R p = ^—\Y B+L+2S , where S is the spin of the par- 
ticle, is multiplicatively conserved. However, from a theoretical point of view, R-parity 
conservation is not needed. Allowing its violation leads to a more general superpotential 
W which can include the following renormalizable gauge invariant additional terms: 

Wal^o — KjkLiLjE k + X' ijk LiQjD k + €iH u Li 
WWo = KjkUiDjDk. 

Here L (Q) are the lepton (quark) doublet superfields, E (U , D) are the lepton (up 
and down quark) singlet superfields, H u is the Higgs superdoublet coupling to up-type 
quarks and leptons, i, j, k are generation indices; q are parameters with dimensions 
while Xijk, A^- fc , X'^ k are dimensionless Yukawa-like couplings. The A^ (A" jfc ) couplings 
are anti-symmetric in the first (last) two indices because of gauge invariance. 

Nevertheless, it is necessary to require that AL ^ and AB ^ terms are not both 
present to avoid a fast proton decay. In this paper, it will be assumed that only one A^ 
coupling is non- vanishing. 

R-parity violation has two major consequences. It allows the decay of the lightest 
supersymmetric particle (LSP), thus discarding it as a candidate to cold dark matter. 
It also allows the supersymmetric particles to be singly produced, via A^ couplings in 
the case of e + e~ collisions. It is this possibility that is explored in this paper. If Am or 
A131 is non-vanishing, a muon sneutrino or a tau sneutrino (an electron sneutrino cannot 
be produced because Am = 0), with spin 0, can be produced in the s-channel [3] (see 
Figure 1). The simplest expression for the cross- section is [4]: 

a(e^e — > v — > X)(s) — 



mi ( s - M?y + M?r? 



A 2 - 

where T(ee) = T(i)j — > e + e~) = —^-M u , j = 2, 3 and T(X) denotes the partial width 

1 (iv; 

for Vj decay to a final state X, with X = e + e~ (sneutrino direct decay), or x^F 
(sneutrino indirect decays). In the indirect decay mode, T(X) is independent of X^. The 
additional t-channel contributions and the interference terms, not shown in this formula 
for reasons of compactness, were included in the signal simulation [5] and final estimation 
of expected events. This cross-section is expected to be very high for M/> ~ (of the 
order of 50 pb for instance for Ayi = 1CT 2 ) and to remain large for masses below the 
centre-of-mass energy, due to initial state radiation effects, again not displayed by this 
formula but taken into account in the analysis. 

Given the present (indirect) upper limits on A m and A131 (Ai 2 i < 0.04 x 100 G e c v/c 2 ' 
A131 < 0.05 x 10 Q j ^ 5 y /e 2 & t 68% confidence level (CL) [6]), the e + e~ decay channel, having a 
cross-section proportional to A 4 , is suppressed compared to the other two (a oc A 2 ), unless 
all neutralinos and charginos are heavier than the sneutrino. The direct decay mode has 
already been investigated by the LEP collaborations [7] by looking for deviations from the 
Standard Model in the cross-sections and asymmetries of e + e~ — > l + l~ . The results were 
presented as an upper limit on Ayi as a function of the sneutrino mass My. The indirect 
decay modes are analysed explicitly here, taking into account any mass and width of the 



2 




Figure 1: Lowest order Feynman diagrams for single sneutrino production and decay: a) 
direct decay, b) indirect decay (vj = u^, u T ; lj — /j,,t; m — 1,4; n = 1,2). 
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sneutrino as a function of the MSSM parameters. The scheme chosen here is a constrained 
MSSM in which the SUSY breaking occurs via gravitational interactions (mSUGRA [1]). 
The relevant parameters are then: M 2 , the SU(2) gaugino mass at the electroweak scale; 
mo, the scalars common mass at the Grand Unified Theories (GUT) scale; /x, the mixing 
mass term of the Higgs doublets at the electroweak scale; and tan /?, the ratio of the 
vacuum expectation values of the two Higgs doublets. The unified trilinear coupling A 
is assumed to be zero. It is also assumed that the running of the A couplings from the 
GUT scale to the electroweak scale does not have a significant effect on the running of 
the gaugino and sfermion masses. 

In this model, the LSP is generally the lightest neutralino Xii with a Am (A131) 
coupling, it decays to e/iz/ e or eev^ (erz/ e or eeu T ). All other sparticles, like the sneutrino, 
can have both direct and indirect decays. In particular, the lightest chargino xt can 
decay either to ee/i, ev e v^ (Am) or to the R p conserving channel XiW*. The latter is 
generally clearly dominant over the former, unless the A coupling is very large or the 
mass difference between the chargino and the neutralino is very small. The decays of the 
heavier neutralinos and charginos are similar, with the additional possibility of longer 
cascade decays always leading to the LSP decay. 

A complete review of the possible decays of the neutralinos and charginos showed 
that the final states of all indirect decays of the sneutrino could be classified into three 
topologies (or channels): 

1) events with two leptons and missing energy; 

2) events with four or six leptons (with or without missing energy); 

3) events with at least two isolated leptons and at least two hadronic jets. 

The first topology mostly comes from v — > x\v, Xi ~^ ^l^v, leading to a 21 + v final 
state. The second one is dominated by v — > xtl~i Xi ~^ Xi^ +z/ > Xi ~^ e ± l T v, leading to 
4/ + 2v. The third topology, being relatively general, often has a high branching ratio. It 
can arise for example from v — > Xi^ \ Xi ~^ XiQQ'i Xi ~^ e ± l T u, leading to 3/ + 2jets+z/, 
or from v — > x\ v -, X2 X1QQ1 Xi ~^ e ± / T z/, leading to 21 + 2jets+2z/. 

The paper is organised as follows. Section 2 lists the data samples that were used in 
the present search for resonant sneutrino production. In section 3, the selection criteria 
are described, as well as the results of the selection. Finally, in section 4 the limits on 
the Aiji couplings are derived by comparing the Standard Model expectations and the 
experimental results. 

2 Data samples 

A detailed description of the DELPHI detector can be found in [8]. The present 
analysis was mainly based on the capability of reconstructing charged particle tracks 
using the tracking devices, particularly the Time Projection Chamber (TPC) but also 
the silicon Vertex Detector, the drift chambers called Inner and Outer Detectors, and the 
forward detectors. The complete system was inside a solenoidal magnetic field of 1.2 T, 
parallel to the beam axis. The analysis also used the lepton identification capabilities 
of the electromagnetic calorimeters (the barrel High density Projection Chamber HPC 
and the Forward Electro-Magnetic Calorimeter FEMC) for the electrons and of the muon 
chambers for the muons. The hadron calorimeter was used to detect neutral hadrons. 

The data from 1997 to 2000 LEP runs were taken at centre-of-mass energies between 
183 GeV and 208 GeV. The registered integrated luminosities after requiring the TPC 
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< v /s>(GeV) 


JC (pb- 1 ) 


182.7 


52.2 


188.6 


153.8 


191.6 


25.1 


195.5 


75.9 


1 on c 


nn o 
OZ.O 


201.6 


43.2 


203.7 


6.3 


205.0 


67.2 


206.5 


78.2 


208.0 


7.3 


All 


592.0 



Table 1: Average centre-of-mass energies and integrated luminosities. 

and all the calorimeters (HPC, FEMC and hadron calorimeter) to be operational are 
given in Table 1. This quality requirement rejected 2% of the luminosity recorded in 
1997 to 1999 and 29% of 2000 integrated luminosity, because one of the twelve sectors 
of the TPC was off before the end of the data taking. Runs with partly inefficient muon 
chambers were kept, and the simulation was adjusted to reproduce the effective efficiency. 

To evaluate the background contamination, different contributions coming from the 
Standard Model were considered. The Standard Model events were produced by the 
following generators: BDKRC [9] for the e + e~ — > e + e~l + l~ four-fermion events of type 77, 
and WPHACT [10] for the other four-fermion events; KK2F [11] for the two-fermion events 
of type e + e~ — >• 7/(7), with / ^ e, r, BHWIDE [12] for the Bhabha events (/ = e) and 
K0RALZ [13] for / = r events; PYTHIA [14] and WPHACT for the 77 — >• hadrons events. 

Signal events were generated with the SUSYGEN 2.20 generator [5]. Samples of 3500 
to 6000 events were generated for nine MSSM parameter sets (see Table 2). 

Simulated events were produced from the generated samples with the standard 
DELPHI simulation program DELSIM [8] and passed through the same reconstruction 
chain as the data. These events were used to design the event selection. 

A faster simulation programme (SGV) [15] was applied to the same generated samples 
in order to be validated for further use in the limits extraction. 

3 Event selection 

In order to select the three categories of final states defined above, a preselection was 
first applied. In the following, charged particles reconstructed from their trajectories in 
the tracking chambers were accepted only if their momenta were greater than 100 MeV/c 
and less than 1.5 times the beam energy. They also had to have a relative momentum 
error less than 100% and impact parameters of at most 4 cm both in the plane transverse 
to the beam and along the beam direction. On the other hand, a neutral particle is 
assumed to be detected as a cluster of energy deposits not associated to a track. In each 
event, the following were required: 

• at least two charged particles; 
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Parameter set 


m 


n 


M 2 








Bri 


Br 2 


Br 3 


1 


150 


175 


125 


1 


37 


68 


0.39 


0.23 


0.38 


2 


190 


275 


155 


1 


65 


116 


0.42 


0.24 


0.34 


3 


200 


195 


175 


1 


64 


106 


0.41 


0.20 


r\ on 

0.38 


4 


207 


125 


385 


0.2 


87 


104 


0.57 


0.15 


0.28 


5 


207 


-75 


155 


0.5 


71 


96 


0.19 


0.24 


0.57 


6 


207 


-125 


115 


1 


63 


127 


0.10 


0.17 


0.72 


7 


207 


305 


135 


1.5 


57 


105 


0.37 


0.24 


0.39 


8 


207 


-285 


65 


2 


36 


81 


0.17 


0.37 


0.46 


9 


220 


285 


185 


1 


80 


142 


0.52 


0.18 


0.31 



Table 2: Values of the SUSY parameters (in GeV/c 2 ) used in the signal simulation, 
and resulting width, masses and branching ratios (Br^ is the branching ratio of channel 
number i). 

• the sum of all charged particle energies greater than 0.1 x ^/s (and in any case, total 
energy of all particles E tot greater than 20 GeV); 

• the total momentum transverse to the beam greater than 5 GeV/c; 

• the absolute value of total electric charge at most 1 if there are less than 7 charged 
particles; 

• at least one charged particle in the barrel (polar angle between 40° and 140°); 

• the absolute value of the cosine of the polar angle of the missing momentum vector 
below 0.95 (or 0.9 in case of exactly two charged particles); 

• at least one isolated (i.e. with no other track in a 5° half-cone centred on its direction) 
identified lepton (electron or muon) with momentum above 5 GeV/c and with a 
maximum angle of 170° with respect to the nearest track. 

The lepton identification used standard DELPHI tools [8]. The electron identification 
algorithm relied on two types of informations: the energy deposited in the electromagnetic 
calorimeters and the dE/dx measurement in the TPC. The muon identification algorithm 
was based on the association of signals in the muon chambers with extrapolated tracks; 
the most efficient set of criteria, as described in [8], was chosen for this analysis. 

The above criteria define the preselection. Then four more series of requirements were 
designed, in order to select the different kinds of topologies. For the two acoplanar 1 
lepton and four or six lepton topologies, the criteria were the same for both X 12 i and Ai 3 i 
couplings; for the semi-leptonic topology, two slightly different selections were applied. 

• For the two acoplanar lepton final states, it was required that: 

— there be exactly two charged particles; 

— not both identified as muons; 

— the acoplanarity be above 40°; 

— the acollinearity 2 be above 50°; 

— the invariant mass of the two leptons 3 be lower than 0.25 x y/s; 

— the angle between the two leptons be lower than 100°. 

The acoplanarity is 180° minus the angle between the transverse momenta of the two charged particles, or of the 
two reconstructed jets (forcing the number of jets to be two with the LUCLUS algorithm [16]) if the event contains more 
particles. 

2 The acollinearity is 180° minus the angle between the two charged particles (or the two jets, see the definition of 
acoplanarity) . 

3 When the second lepton was not identified, it was assumed to be an electron. 
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• For the four or six lepton final states, the requirements were: 

— four or six charged particles; 

— at least two identified leptons (electrons or muons); 

— the resolution parameter of the Durham algorithm [17] for which the event 
changes from four to three jets 2/34 greater than 10~ 4 . 

• For the semi-leptonic final states, A m coupling: 

— at least 7 charged particles and at most 25; 

— at least two identified leptons (electrons or muons); 

— the transverse momentum of the second most energetic lepton had to be above 
0.05 x v/i; 

— y 3 4 had to be greater than 10 -3 ; 

— when the number of jets was forced to four, at least two of them were required 
to be thin, that is to have a total (track + neutral) multiplicity not exceeding 4. 

• For the semi-leptonic final states, A131 coupling: 

— at least 7 charged particles and at most 25; 

— at least two identified leptons (electrons or muons), including at least one iden- 
tified electron; 

— ?/34 had to be greater than 10 -3 ; 

— when the number of jets was forced to four, at least two of them were required 
to be thin; 

— the missing energy (y/s — E tot ) had to be greater than 0.25 x y/s. 

The numbers of data and of SM Monte Carlo events after the preselection are shown 
in Table 3. Distributions of some important variables are shown in Figures 2 and 3 at 
the preselection level. The plots were chosen so as to give at least an example of each 
centre-of-mass energy. The agreement between real data and simulated SM background 
is good. 

Some examples of signal distributions are also given, for the following parameters: 
X in = 0.05, m = y/s, tan/3 = 1.5, /1 = -125 GeV/c 2 , M 2 = 115 GeV/c 2 . The missing 
energy signal is given for the semi-leptonic channel, j = 3, and scaled by a factor of 5; 
the second lepton transverse momentum and the thin jet multiplicity signals are given 
for the semi-leptonic channel, j = 2, and scaled by a factor of 5 and 2 respectively; the 
log 10 2/ 34 signal is given for the four-lepton channel, j = 3, and scaled by a factor of 20; 
the acollinearity and the two-lepton angle signals are given for the two-lepton channel, 
j = 2, and scaled by a factor of 20; the acoplanarity and the two-lepton invariant mass 
signals are given for the two-lepton channel, j = 3, and scaled by a factor of 20 and 50 
respectively. 

The efficiency of the selections, including the preselection, depends on the SUSY pa- 
rameters. It is shown in Figure 4 for the nine parameter sets which were fully simulated 
with DELSIM (Table 2) using a A m coupling and y/s = 206.5 GeV. The Figure also shows 
the efficiencies evaluated with SGV; they are compatible. The fact that SGV efficiencies 
were systematically lower than DELSIM efficiencies in the four or six lepton channel was not 
compensated for, first because this has a small impact on the global efficiency and second 
because it gives conservative results. The efficiencies for a A131 coupling are comparable, 
although always lower. 

The expected background at the end of the selections is mainly composed of four- 
fermion events. The 77 background is totally negligible. The background coming from 
the two-fermion events is small and its proportion decreases when the centre-of-mass 
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(GeV) 


SM 


Data 


182.7 


1003±11 


1071 


188.6 


2909±16 


2913 


191.6 


466±3 


484 


i nr r 

195.5 


1387±8 


1466 


199.5 


1487±8 


1537 


201.6 


772±5 


777 


203.7 


110±1 


99 


205.0 


1180±7 


1210 


206.5 


1345±6 


1322 


208.0 


126±1 


115 



Table 3: Number of events after the preselection (SM= Monte Carlo simulated SM back- 
ground normalised to the data luminosity). Errors are statistical only. 

energy increases. For instance, it is 9% in the two-lepton channel at 183 GeV and 4% at 
206.5 GeV. 

The number of data and of SM Monte Carlo events at the end of the selections is 
shown in Table 4. There is no significant excess of data in any of the three channels and 
in any of the centre-of-mass energy samples. 

4 Limits on Am and A131 couplings 

Besides being used as an event generator, SUSYGEN was also used to scan a wide part 
of the MSSM parameter space and compute all the cross-sections of the signal, with 
Ayi = 5 x 10~ 3 . In the model adopted for this search and described in section 1, all 
SUSY phenomenology can be derived from the four parameters tan/5, mo, M2 and /x, 
plus the centre-of-mass energy for the kinematics. The parameter sets explored in the 
scan were: 

• y/s= 182.7, 188.6, 199.5 and 206.5 GeV, 

• tan/? = 1.5 or 30, 

• m = 100 to 230 GeV/c 2 (170 to 215 GeV/c 2 in steps of 1 GeV/c 2 , 100 to 170 GeV/c 2 
in steps of 10 GeV/c 2 , 215 to 230 GeV/c 2 in steps of 5 GeV/c 2 ), 

• M 2 = 5 to 405 GeV/c 2 in steps of 10 GeV/c 2 , 

• n = -305 to 305 GeV/c 2 in steps of 10 GeV/c 2 . 

The cross-sections at yfs = 191.6 GeV were taken from the y/s = 188.6 GeV scan 
simply assuming y/s + 3 GeV; in the same way, the cross-sections at 195.5 and 201.6 GeV 
were taken from the 199.5 GeV scan and the cross-sections at 203.7, 205.0 and 208.0 GeV 
were taken from the 206.5 GeV scan assuming the corresponding centre-of-mass energy 
shifts. 

The sneutrino mass Mc was assigned the value of m , thus slightly departing from a 
strict mSUGRA model. This is a conservative hypothesis since the sneutrino mass tends 
to decrease whereas the gaugino masses are left untouched. 

Very small values of the sneutrino total width (r^ < 150 MeV/c 2 ) correspond to 
regions of the parameter space where the sneutrino is lighter than the gauginos; they can 
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Figure 2: Examples of data-simulation comparison at preselection level. The dashed lines 
are examples of signal distributions (see text). 
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Figure 3: Additional examples of data-simulation comparison at preselection level (the 
log 10 2/34 and jet multiplicity distributions are given for events with at least four charged 
particles). The dashed lines are examples of signal distributions (see text). 
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Figure 4: SGV-DELSIM comparison for nine SUSY parameter sets with different sneutrino 
masses and total widths (see Table 2) in the case of a A m coupling and for tan/3 = 1.5, 
y/s = 206.5 GeV. For this comparison, the same number of events was simulated with 
SGV and DELSIM. 
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Table 4: Number of events at the end of the selection (SM = Monte Carlo simulated 
SM background). Errors are statistical only. Note that the sum of all channels is not 
independent for A m and A131 due to the common analysis for the first two channels. 

hardly be detected by the present analysis, however they are covered by e + e~ — > l + l~ 
analyses [7]. 

Regions of the parameter space already excluded by the precision measurements 
at LEP1 were not further explored in the present scans. This condition was imple- 
mented in the following way. Using the expression for the cross-section at the resonance, 
a(e+e- -> Z -> X) = 12tt^|, the limit T new < 6.6 MeV/c 2 at 95% CL [18] can be 

converted into an upper limit on the cross-section of new decay modes: a new < 157.2 pb 
at 95% CL. SUSY parameter sets for which the total cross-section of pair production 
of charginos and neutralinos at \/s ~ M z was larger than this limit were considered as 
excluded by LEP1. 

To derive the limits on one A coupling, each centre-of-mass energy was first considered 
separately. The three channels being totally independent due to the charged particle 
multiplicity criterion, they were summed up (Table 4). On the other hand, the (M/>,rV) 
plane was divided into rectangular bins of size (1 GeV/c 2 , 50 MeV/c 2 ). For each set of 
parameters entering a given bin, the output of two scans was used. First the SUSYGEN 
scan to get a, the total cross-section expected for e + e~ — > v — > X, X representing all 
final states from indirect sneutrino decays. From the expression given in the introduction, 
a = o" x A 2 . cr is almost independent of A as long as A is reasonably small: for A = 0.1, 
do is at most a few percent lower than for A = 10" 4 . Second a SGV scan on the same 
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parameters (however with wider steps, of 10 GeV on the full m Q range and of 20 GeV for 
Li and M 2 ), to obtain the global efficiency of the analysis, e, combining the three separate 
channel efficiencies according to the branching ratios predicted by SUSYGEN. 1000 events 
were simulated for each parameter set. 

An upper limit at 95% CL on the total number of signal events (N up ) compatible with 
the data and expected background was then calculated for each (M^, r„) bin, combining 
the ten centre-of-mass energy samples considered as independent samples. This was done 
using the Bayesian method described in reference [19]; the relative probabilities of each 
centre-of-mass energy were taken as 



Wi 



10 ' 
3=1 

where (<T e)i is the lowest value of such a product in the considered bin and J Ci is the 
integrated luminosity of the i th sample. The 95% CL upper limit on A in each bin was 
then 

A < 



TV 



10 



\ 3=1 



The whole procedure was repeated for the second coupling. 

The results are shown in Figures 5 to 8 corresponding to the two values of tan/3 
considered. In a very large fraction of the (Mp,Fp) plane the obtained upper limit on 
Aiji is at most 0.1, and below 0.01 in still the major part of the area allowed by the 
parameters range. 

There are two main sources of systematic errors on these results. One is the es- 
timation of the expected SM background. Uncertainties arise from the Monte Carlo 
statistics (at most ±5%, see Table 4), from the detector response simulation and from 
the cross-sections evaluation and event modelling. In the last case, a comparison of dif- 
ferent generators gave at most ±5% difference, mainly coming from the hadronisation 
modelling, for the four-fermion processes which are the dominant background. This un- 
certainty is smaller for two-fermion processes and larger for 77 events, which however do 
not affect the final results because very few events of that kind remain at the end of the 
selection. The resulting systematic uncertainty on the limit on A, evaluated separately 
for each centre-of-mass energy, is of the order of ±2%. The second source of system- 
atic errors is the estimation of the selection efficiency, performed with a fast simulation 
and however statistically limited (up to ±5% on low branching ratio channels). Here an 
important uncertainty arises from the track reconstruction and the lepton identification 
efficiencies. Discrepancies between DELSIM and SGV are at the level of ±4% for the global 
efficiency (Figure 4). The resulting variation of the limit on A, using the prescription of 
reference [20], is ±3% for the sample with highest luminosity, and at most ±1% for the 
other samples. These effects were considered small and were not included in the plots of 
Figures 5 to 8. 

In order to make the results easier to read, the upper limits on the A couplings were 
also derived as a function of the sneutrino mass only. This was simply done by keeping 
the most conservative limit for each sneutrino mass, assuming > 150 MeV/c 2 . The 
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results can be seen in the same Figures as for the two-dimensional limits. The limits 
clearly show the centre-of-mass energy structure of the data samples; they are especially 
stringent for My ~ yfs. 

5 Conclusion 

The possibility of single production of supersymmetric particles was explored but none 
was seen. Upper limits at the 95% CL were derived in the mSUGRA constrained MSSM 
framework with low and high values of tan (3 for the two possible R p violating couplings 
A121 and Ai 3 i. They are at the level of 2 to 3xl0 -3 , depending on Mo, when it is close 
to the centre-of-mass energy. They are slightly better for high tan (3 and slightly worse 
for A131 as compared to A^i- In any case, they are one or two orders of magnitude better 
than the published indirect limit, even for sneutrino masses outside the centre-of-mass 
energies of the analysed data, due to the i-channel and, more importantly, to the initial 
state radiation below the resonance. 

Acknowledgements 

We are greatly indebted to our technical collaborators, to the members of the CERN- 
SL Division for the excellent performance of the LEP collider, and to the funding agencies 
for their support in building and operating the DELPHI detector. 
We acknowledge in particular the support of 

Austrian Federal Ministry of Education, Science and Culture, GZ 616.364/2-III/2a/98, 
FNRS-FWO, Flanders Institute to encourage scientific and technological research in the 
industry (IWT), Belgium, 

FINEP, CNPq, CAPES, FUJB and FAPERJ, Brazil, 
Czech Ministry of Industry and Trade, GA CR 202/99/1362, 
Commission of the European Communities (DG XII), 
Direction des Sciences de la Matiere, CEA, France, 

Bundesministerium fur Bildung, Wissenschaft, Forschung und Technologie, Germany, 
General Secretariat for Research and Technology, Greece, 

National Science Foundation (NWO) and Foundation for Research on Matter (FOM), 
The Netherlands, 
Norwegian Research Council, 

State Committee for Scientific Research, Poland, SPUB-M/CERN/PO3/DZ296/2000, 

SPUB-M/CERN/PO3/DZ297/2000, 2P03B 104 19 and 2P03B 69 23(2002-2004) 

JNICT- Junta Nacional de Investigagao Cientifica e Tecnologica, Portugal, 

Vedecka grantova agentura MS SR, Slovakia, Nr. 95/5195/134, 

Ministry of Science and Technology of the Republic of Slovenia, 

CICYT, Spain, AEN99-0950 and AEN99-0761, 

The Swedish Natural Science Research Council, 

Particle Physics and Astronomy Research Council, UK, 

Department of Energy, USA, DE-FG02-01ER41155. 



14 



References 

[1] H.P. Nilles, Phys. Rep. 110 (1984) 1; 

H.E. Haber and G.L. Kane, Phys. Rep. 117 (1985) 75. 
[2] P. Fayet, Phys. Lett. B69 (1977) 489; 

G. Farrar and P. Fayet, Phys. Lett. B76 (1978) 575. 
[3] S. Dimopoulos and L.J. Hall, Phys. Lett. B207 (1988) 210. 
[4] V. Barger et al., Phys. Rev. D40 (1989) 2987. 

[5] S. Katsanevas and P. Morawitz, Comput. Phys. Commun. 112 (1998) 227. 

[6] H. Dreiner, in Perspectives on Super symmetry, edited by G.L. Kane (World Scientific, 

Singapore, 1998) p. 462. 
[7] R. Barate et al, ALEPH coll., Eur. Phys. J. C12 (2000) 183; 

P. Abreu et al., DELPHI coll., Phys. Lett. B485 (2000) 45; 

M. Acciarri et al, L3 coll., Phys. Lett. B489 (2000) 81; 

G. Abbiendi et al, OPAL coll., Eur. Phys. J. C13 (2000) 553. 
[8] P. Abreu et al., DELPHI coll., Nucl. Inst, and Meth. A378 (1996) 57; 

P. Aarnio et al, DELPHI coll., Nucl. Inst, and Meth. A303 (1991) 233. 
[9] F.A. Berends, P.H. Daverveldt and R. Kleiss, Comput. Phys. Commun. 40 (1986) 

271. 

[10] E. Accomando and A. Ballestrero, Comput. Phys. Commun. 99 (1997) 270; 

E. Accomando, A. Ballestrero and E. Maina, hep-ph/0204052, to appear in Comput. 
Phys. Commun. 

[11] S. Jadach, B.F.L. Ward and Z. Was, Comput. Phys. Commun. 130 (2000) 260. 

[12] S. Jadach, W. Placzek and B.F.L. Ward, Phys. Lett. B390 (1997) 298. 

[13] S. Jadach, B.F.L. Ward and Z. Was, Comput. Phys. Commun. 79 (1994) 503. 

[14] T. Sjostrand, Comput. Phys. Commun. 39 (1986) 347. 

[15] R. Keranen et al, Eur. Phys. J. direct C7 (2000) 1. 

[16] T. SjSstrand, Comput. Phys. Commun. 82 (1994) 74. 

[17] S. Catani et al., Phys. Lett. B269 (1991) 432. 

[18] P. Abreu et al, DELPHI coll., Eur. Phys. J. C16 (2000) 371. 

[19] V.F. Obraztsov, Nucl. Inst, and Meth. A316 (1992) 388 and erratum A399 (1997) 
500. 

[20] R.D. Cousins and V.L. Highland, Nucl. Inst, and Meth. A320 (1992) 331. 



15 




120 140 160 180 200 220 

Sneutrino mass(GeV/c 2 ) 




120 140 160 180 200 220 



Sneutrino nnass(GeV/c 2 ) 

Figure 5: For tan/3 = 1.5, upper limit on A i2 i as a function of M„ and Y„ (top) and 
as a function of assuming Yy > 150 MeV/c 2 (bottom). The white zone in the top 
plot corresponds to non existing sneutrino widths given the /i parameter range. The area 
entitled 'no limit on A' corresponds to upper limits larger than 0.1. The indirect limit 
coming from precision measurements is drawn in the bottom plot assuming Mg R = M„. 
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Figure 6: For tan/3 = 30, upper limit on Ai 2 i as a function of M 5 and r c (top) and as 
a function of M„ assuming T„ > 150 MeV/c 2 (bottom). The white zone in the top plot 
corresponds to non existing sneutrino widths given the \x parameter range. The area 
entitled 'no limit on A' corresponds to upper limits larger than 0.1. The indirect limit 
coming from precision measurements is drawn in the bottom plot assuming Mg R = M„. 
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Figure 7: For tan/3 = 1.5, upper limit on A131 as a function of M„ and Y„ (top) and 
as a function of assuming Yy > 150 MeV/c 2 (bottom). The white zone in the top 
plot corresponds to non existing sneutrino widths given the /i parameter range. The area 
entitled 'no limit on A' corresponds to upper limits larger than 0.1. The indirect limit 
coming from precision measurements is drawn in the bottom plot assuming Mg R = M„. 
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Figure 8: For tan/3 = 30, upper limit on A131 as a function of M 5 and r c (top) and as 
a function of M„ assuming T„ > 150 MeV/c 2 (bottom). The white zone in the top plot 
corresponds to non existing sneutrino widths given the \x parameter range. The area 
entitled 'no limit on A' corresponds to upper limits larger than 0.1. The indirect limit 
coming from precision measurements is drawn in the bottom plot assuming Mg R = M„. 



